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ABSTRACT 


This  thesis  investigates  the  problem  of  economic 
scheduling  of  an  interconnected  hydro-thermal  power 
system,  over  short  intervals.  The  minimization  of  the 
operating  cost  of  the  system  is  treated  as  a  Lagrange 
problem  using  the  Calculus  of  Variations.  A  new  method 
has  been  developed  for  obtaining  scheduling  equations 
for  the  hydro-plants  of  the  system.  The  effect  of  the 
variation  of  head  upon  the  hydro-plant  characteristics 
and  the  transmission  line  losses  are  taken  into  consi¬ 
deration. 

The  scheduling  equations  for  the  hydro-plants 
located  on  separate  streams  are  shown  to  be  equivalent 
to  those  derived  by  Glimn  and  Kirchmayer  [14] . 

New  scheduling  equations  are  developed  in  this 
thesis  for  the  common-flow  hydro-plants  of  the  system. 

For  simplicity  it  has  been  assumed  that  only  two  hydro¬ 
plants  are  located  on  the  same  stream.  The  time  taken 
by  water  to  flow  from  the  upstream  plant  to  the  down¬ 
stream  plant  is  taken  into  account.  The  method  is  general 
in  nature  and  can  be  extended  to  cases  where  more  than 
two  hydro-plants  are  located  on  the  same  stream. 

The  scheduling  equations  are  applied  to  a  simple 
power  system  to  test  their  effectiveness.  The  computing 
technique  used  for  doing  so  is  also  discussed. 
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CHAPTER  I 


INTRODUCTION 


1-1  General 

The  generation  of  electrical  power  by  thermal  and 
hydraulic  means,  using  the  natural  resources  such  as  water 
and  coal,  represents  a  very  important  factor  in  a  nation's 
economy.  The  effective  utilization  of  these  natural  resour¬ 
ces,  for  the  generation  of  electrical  power,  will  reduce 
not  only  the  operating  costs,  but  also  future  capital 
commitments  for  additional  generating  equipment.  Relatively 
small  deviations  from  the  optimal  coordination  of  the 
expenditure  of  these  resources  can  cause  considerable 
monetary  losses  over  a  period  of  time. 

Thus,  the  problem  is  basically  one,  in  which  the 
limited  available  natural  resources  have  to  be  most  effec¬ 
tively  utilized  to  meet  certain  load  requirements.  The 
generation  schedules  for  each  of  the  thermal  and  hydro¬ 
electric  plants  should  be  such  that  the  cost  of  operating 
the  system,  over  the  given  optimization  period,  is  a 
minimum. 

The  methods  of  scheduling  the  operation  of  pure 
thermal  systems  [1,2,3],  by  using  the  incremental  cost 
techniques,  have  been  widely  discussed  in  the  literature 
and  have  been  put  to  extensive  use.  The  problem  of 
preparing  an  economic  operating  schedule  for  a  hydro- 
thermal  system  is  more  complex.  The  essential  difference 
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between  the  operation  of  an  all  thermal  and  a  hydro- 
thermal  system  is  that  the  operation  of  the  thermal 
plants  depends  only  upon  the  conditions  which  exist 
from  instant  to  instant,  while  the  operation  of  the 
hydrothermal  system  depends  upon  the  conditions  which 
exist  over  the  entire  optimization  interval  [4]. 

The  hydro-thermal  systems  can  be  divided  into  two 
catagories.  In  the  first,  the  hydraulic  plants  are 
located  on  separate  streams,  and  in  the  second  some  or 
all  the  hydraulic  plants  are  situated  on  the  same  stream. 
In  the  latter  case,  the  discharges  from  the  upstream 
plants  affects  the  operation  of  the  downstream  plants 
as  well  [5].  In  this  case,  the  time  taken  by  the  water 
to  flow  from  the  upstream  plants  to  the  downstream  plants 
should  also  be  taken  into  consideration. 

The  optimization  problem  can  be  classified  on  the 
basis  of  the  duration  of  the  optimization  period: 

(a)  Long-range  problem  -  The  period  of  optimization 
may  range  from  a  few  months  to  one  year. 

(b)  Short-range  problem  -  The  optimization  interval 
may  extend  from  24  hours  to  a  week. 

The  difficulties  in  predicting  the  natural  inflows, 
makes  the  long  range  problem  rather  complex.  As  a  matter 
of  fact,  only  probabilistic  techniques  may  give  worthwhile 
results  for  a  long-range  problem.  It  is  relatively  easier 
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to  predict  the  natural  inflows  over  shorter  intervals. 

A  technique  similar  to  the  Basic  Rule  Curve  [6] ,  may  be 
used  to  decide  the  amount  of  water,  which  is  to  be  used 
over  the  optimization  interval  of  24  hours.  The  short- 
range  problem  is  the  subject  of  investigation  in  this 
thesis . 

The  short-range  problem  has  been  attacked  by 
other  investigators  in  several  ways.  In  almost,  all 
cases  a  mathematical  model  of  the  system  is  constructed 
and  solved.  The  method  of  formulating  the  problem 
depends  upon  the  mathematical  technique  adopted  for 
solving  it.  The  Variational  Calculus  [7,8,9],  the 
Maximum  Principle  [9] ,  Dynamic  Programming  [9]  and  the 
Gradient  Methods  [9]  have  been  successfully  used.  Dynamic 
Programming  and  the  Gradient  Methods,  for  a  given  set 
of  conditions,  actually  compute  all  possible  combinations 
of  the  thermal  and  hydro-plant  generations  from  hour  to 
hour;  and  then  select  that  combination  which  results  in 
minimum  cost  over  the  optimization  interval. 

The  methods  of  Calculus  of  Variations,  and  the 
Maximum  Principle  depend  upon  the  iterative  calculations 
which  converge  to  the  desired  answers.  The  feasibility 
of  using  Pontryagin's  Maximum  Principle  for  preparing 
the  optimal  generation  schedules,  has  been  shown  by 
Dhalin  and  Shen  [10]  ,  and  Hano,  Tamure  and  Narita  [11] . 
But  as  the  number  of  the  plants  in  the  system  increases, 
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it  becomes  quite  difficult  to  match  the  two  point  boundary 
conditions  for  the  state  and  the  adjoint  equations  and  to 
evaluate  the  Hamiltonian  numerically. 

The  Variational  techniques  have  been  extensively 
discussed  and  used  to  obtain  the  optimal  generation  schedules. 
Ricard  [12]  used  the  Variational  Calculus  to  arrive  at  the 
well-known  Ricard's  equation.  Cypser  [5]  posed  the  problem 
using  the  Calculus  of  Variations  and  solved  the  resulting 
Euler  equations  by  the  Steepest  Descent  method.  He  discused 
the  common-flow  problem  as  well.  Chandler,  Dandeno,  Glimn  and 
Kirchmayer  [13]  advanced  the  concept  of  incremental  cost  to 
combiined  hydro-thermal  systems,  by  treating  the  hydro-plants 
as  constant  head  plants.  Glimn  and  Kirchmayer  [14]  improved 
upon  this  method  by  considering  the  effect  of  variation  of 
head  upon  the  hydro-plant  characteristics.  Arismunandar  [1] 
made  an  effort  to  formulate  the  problem  for  a  combined  hydro- 
thermal  electric  system,  but  made  no  effort  to  consider 
systems  containing  the  Common-flow  plants.  Menon  [15] 
used  the  Calculus  of  Variations  approach  to  formulate  the 
problem  and  then  used  the  Euler  equations  for  constructing 
sets  of  minimizing  sequences. 
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1-2  Scope  of  the  Thesis 

The  purpose  of  this  investigation  is  to  obtain, 
using  Calculus  of  Variations,  general  scheduling  equations 
for  a  hydro-thermal  power  system  containing  common- 
flow  hydro-electric  plants.  The  integral  of  the  cost  of 
operating  the  thermal  plants  of  the  system,  over  the 
optimization  interval,  is  taken  as  the  performance  index 
for  the  problem.  The  transmission  losses  in  the  system 
and  the  effect  of  variation  of  head  upon  the  hydro-plant 
characteristics  are  taken  into  consideration. 

The  scheduling  equations  for  the  hydro-plants  located 
on  separate  streams,  are  shown  to  be  ecruivalent  to  those  by 
Glimn  and  Kirchmayer  [14]  (their  equations  are  extensions 
of  Ricard's  equations).  New  coordinating  equations  are 
obtained  for  plants  located  on  the  same  stream. 

An  example  to  illustrate  this  method  is  discussed. 

The  system,  under  study,  consists  of  one  thermal  and  two 
hydro-plants.  The  hydro-plants  have  storage  facilities 
on  the  same  river.  The  transmission  losses  and  the  time 
lag  between  the  reservoirs  are  considered  negligible.  The 
scheduling  equations  for  the  common-flow  plants  are  simpli¬ 
fied  for  the  given  system. 

The  numerical  results  are  obtained  with  the  help 
of  the  IBM-360  digital  computer  system.  Fortran  IV 
programming  language  is  used.  The  programme  is  discussed 


in  some  detail. 
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CHAPTER  II 

SYSTEM  CHARACTERISTICS  AND  CONSTRAINTS 

2-1  A  General  Statement  of  the  Problem 

Consider  a  system  containing  M  thermal  and  N 
hydro-electric  plants,  interconnected  electrically 
by  a  network  of  transmission  lines.  A  predetermined 
limited  amount  of  water  is  to  be  used  by  each  hydro¬ 
plant,  over  the  optimization  period.  The  natural  inflows 
and  the  load  demand  are  known  functions  of  time.  It  is 
desired  to  operate  this  system  to  meet  the  load  demand  in 
such  a  way  that  the  total  operating  cost  of  the  system, 
over  the  optimization  interval,  is  a  minimum.  Some  or 
all  the  hydro-plants  may  be  located  on  the  same  stream. 

The  transmission  losses  and  the  effect  of  variation  of 
head  upon  the  hydro-plant  characteristics  are  to  be  taken 
into  consideration. 

In  order  to  construct  a  suitable  mathematical  model 
of  the  system,  the  characteristics  of  the  various  elements 
of  the  problem  should  be  examined  in  detail. 

The  main  elements  of  the  problem  are: 

(i)  Hydro-electric  power  plants, 

(ii)  Thermal  plants  and  the  cost  function, 

(iii)  Transmission  losses  in  the  system,  and 

(iv)  The  system  restrictions. 
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2-1  Hydro-electric  Power  Plants 


The  presence  of  the  hydraulic  plants  in  a  system 
makes  the  problem  quite  complex.  Inaccuracies  in  predic¬ 
ting  the  natural  inflows  to  the  reservoirs,  further  compli¬ 
cates  the  problem.  In  the  following  analysis  it  is  assumed 
that  the  natural  inflows  to  the  reservoirs  of  the  hydro¬ 
electric  plants  are  known  functions  of  time. 

It  is,  however,  interesting  to  compare  the  perfor¬ 
mance  of  the  steam  and  hydraulic  power  plants.  A  hydraulic 
turbine  can  go  from  a  stationary  position  to  full  load  in 
1-3  minutes  [16],  while  it  may  take  hours  to  fire  up  a 
steam  plant  and  put  it  on  line.  A  hydraulic  turbine 
spinning  under  no  load  will  pick  up  load  almost  instanta¬ 
neously.  The  hydraulic  plants  are  inexpensive  means  of 
reliable  reserves  and  are  quite  adaptive  to  load  variations. 

The  output  of  a  hydraulic  turbine  is  given  by: 


P 


H 


q.h.  n 


Kw 


11.8 


(2.1) 


where 

P„  =  Output  of  the  turbine  in  Kw 
ti 

3 

q  =  discharge  in  ft  /sec 

h  =  net  head  in  ft 

n  =  efficiency  of  the  turbine,  under  the 
operating  conditions. 
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The  efficiency  (n)  of  a  turbine  depends  upon  its 

design  and  the  operating  condition.  Thus,  for  a  given 

turbine,  n  can  be  represented  as  a  nonlinear  function  of 

the  discharge  and  the  net  head.  A  hydro-electric  power 

plant  consists  of  a  number  of  such  turbines.  In  general, 

the  output  of  a  power  plant  can  be  represented  as  a 

function  of  the  head  and  the  discharge.  This  is  done  by 

collecting  the  input-output  data  for  a  plant  and  fitting 

surfaces  into  this  experimental  data  [14].  Thus,  the 

th 

output  of  the  j  hydroplant  (at  any  instant  of  time  t) 
can  be  written  as: 

PHj  =  PHj  V  tJ  <2‘2) 

The  problem  stated  in  section  2-1,  is  posed  as 
a  fixed  end  point  Lagrange  problem  in  the  Calculus  of 
Variations  in  the  subsequent  chapters.  To  do  this,  it  is 
desirable  to  represent  the  output  of  the  hydro-electric 
plants  as  a  function  of  storage  at  the  plant  and  the  rate 
of  change  of  storage  at  the  plant  and  at  all  the  upstream 
plants  (if  any)  [5].  This  is  done  by  representing  the 
discharge  and  the  net  head  at  the  plant  as  functions  of 
the  storages  and  their  time  differentials  as  explained 
in  the  following  paragraphs. 
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In  the  following  analysis,  it  is  assumed  that  only 
the  k^and  k+lst  plants  are  located  on  the  same  stream 


(both  plants  have  adequate  reservoir  capacity) ,  the  rest 
of  the  plants  are  located  on  separate  streams. 

(a)  Discharge 

(i)  Plants  on  separate  rivers:  The  discharge  can 
be  represented  as  a  function  of  the  natural 
inflow  to  the  reservoir  of  the  plant  and  time 
differential  of  the  storage  at  the  plant.  Thus 


for  the  jth  plant 


(2.3) 


F, (t)  -  s . (t) 

J  J 


(ii)  Common- flow  plants:  When  the  hydroplants  are 
located  on  the  same  stream,  the  discharges  of 
the  upstream  plants  affect  the  performance  of 
the  downstream  plants  as  well.  The  water 


th 

released  from  the  upstream  (k  hydro-plant ) plant 


takes  time  t  to  reach  the  downstream  plant 
s  t 

(k+1  hydroplant);  that  is,  the  water  released 

th 

at  an  instant  (t-x)  from  the  k  plant  shall 


reach  the  k+lst  hydroplant  at  an  instant  t. 

th, 

Hence,  for  the  k  plant 


Fk(t)  '  sk(t) 


(2.4) 


:■  <  :  •.  r.  s 
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s  t 

and  for  the  k+1  plant 

qk+i(t)  =  Fk+i(t)  +  qk(t_T)  ■  sk+i(t)  (2-5) 


The  above  expressions  can  easily  be  extended  to  a 
case  in  which  more  than  two  plants  are  located  on  the 
same  stream. 


(b)  Net  head 

The  available  head  for  power  generation  is  given 

by  the  difference  between  the  forebay  and  the  tailrace 

elevation  minus  the  losses  due  to  the  friction, etc .  in 

th 

the  conduits.  That  is,  net  head  at  the  j  hydroplant 


(2.6) 


The  forebay  elevation  can  be  represented  as  a 
function  of  storage  in  the  reservoir  of  the  plant  while 
the  tailrace  elevation  and  the  frictional  losses  can  be 
represented  as  functions  of  discharge.  The  discharge, 
however,  can  be  represented  as  a  function  of  time 
differentials  of  the  storages  as  discussed  in  the 
previous  paragraphs.  Thus  the  net  head  can  be  represented 
as  a  function  of  storage  at  the  plant  and  of  the  time 
differential  of  the  storage  at  the  plant  and  at  the 
upstream  plants  (if  any) . 
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(c)  Output  of  the  Hydro-electric  Plants 

The  output  of  a  hydroplant  is  a  function  of 
discharge  and  net  head  at  the  plant,  as  discussed  while 
deriving  the  equation  (2.2).  The  discharge  and  the 
net  head  can  be  represented  as  a  function  of  the  storage 
and  the  time  differential  of  the  storages  and  the  natural 
inflows.  As  the  natural  inflows  are  uncontrollable  variables, 
and  since  their  values  at  various  instants  are  known,  the 
plant  outputs  can  be  represented  as  functions  of  the 
storages  and  their  time  differentials  only. 


(i)  Plants  on  separate  rivers: 


P 


Hj 


(2.7) 


(ii)  For  the  common- flow  plants: 


The  output  of  the  k  (upstream)  plant  is 


given  by: 


P 


Hk 


PHk  (sk'  5k' 


(2.8) 


s  t 

and,  the  output  of  the  k+1  (downstream) 


plant  is  given  by 


PHk+l(t)  “  PHk+l(sk+l 


(t) ,  (qk (t-x) -sk+1 (t) ) ) 


PHk+l (sk+l (t) 


(sk  (t— x ) +sk+1 (t) ) )  (2.9) 


f  P  •  v. 19b 


' 

■ 
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sk (t-T )  can  be  represented  as  some  nonlinear  function  of 


sk(t)  and  t 


sk (t-T )  =  f (sk (t)  ,  t ) 


(2.10) 


If  t  is  considered  to  be  a  constant,  then 


PHk+l(t)  PHk+l^Sk+l 


(t) f  f (sk  (t) )  ,  sk+1 ( t) ) 


(2.11) 


2-3  Thermal  Plants  and  the  Cost  Functions 

The  problem, under  study,  is  one  in  which  the  cost 
of  operating  the  system  over  a  given  interval  is  to  be 
minimized,  while  it  is  desired  to  meet  a  given  load  demand. 
The  operating  costs  are  only  those  costs  which  are  directly 
related  to  the  output  of  the  generating  units.  The  capital 
costs  on  the  reservoirs  and  the  generating  units  are  not 
taken  into  account  here, since  these  are  not  functions  of 
the  power  outputs .  Thus  the  cost  of  operating  the  system 
is  attributed  to  the  cost  of  operating  the  thermal  plants 
only.  The  operating  cost  [17  ]  of  a  thermal  plant  consists 
of : 

(i)  Starting  cost, 

(ii)  No  load  spinning  cost, 

(iii)  The  loading  cost. 


■ 

. 

■ 

3+  .  •  {  *. ’ 


• 
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The  starting  cost  depends  upon  the  down-time  of  the 
units.  If  a  unit  is  stopped  and  then  restarted  immediately, 
there  will  be  no  or  little  starting  cost.  In  general. 


$ 


starting 


(1-e 


(2.12) 


where 

a  =  cooling  rate 

$Q  =  cold  start  cost 

The  no  load  and  the  loading  cost  indicate  the  fuel 
expenditure  to  run  the  unit  once  it  has  been  started  and 
placed  on  line.  The  no  load  cost  is  a  constant  rate  of 
fuel  expenditure  incurred  as  long  as  the  unit  is  running, 
while  the  loading  costs  are  directly  related  to  the  output 
of  the  thermal  plants. 

By  making  suitable  assumptions,  the  cost  of  operation 
of  a  thermal  plant  can  be  represented  as  a  quadratic  or 
cubic  polynomial  of  the  output  of  the  plant. 

Thus,  for  the  i*'*1  thermal  plant, 

CTi  =  CTi(PTi,  t)  (2.13) 

r  —  1,  2,  ....,  M 

The  cost  function  for  the  whole  of  the  system  can 
be  represented  as  a  sum  of  the  cost  functions  of  all  the 
thermal  plants  in  the  system. 


f  3£/S  'M\J  .  o  '  1  ll  -:c  ,  >ol 

io  >i!  i  t  .itdoo  s 

•  .:  ;  -  'V  <r:  .  , 
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That  is. 


M 

Z 

i=l 


C.  (P 


Ti 


t) 


(2.14) 


or 

Cp  —  (Ppi  /  Pp2  /  Ppj^  /  t )  (2.15) 


2-4  Transmission  Line  Losses 


Due  to  the  development  of  large  integrated  power 
systems  and  the  interconnections  between  networks,  it  is 
necessary  to  consider  not  only  the  cost  of  operating  the 
generating  units  but  also  the  cost  of  transmitting 
electrical  power  from  the  generating  stations  to  the 
loads.  This  is  taken  care  of  by  considering  the  power 
losses  in  the  transmission  lines  while  formulating  the 
problem. 

From  the  point  of  view  of  mathematical  formulation, 
it  is  desirable  to  represent  the  power  losses  in  the 
system  as  a  function  of  the  source  outputs.  It  is  done 
by  using  the  approximate  transmission  loss  formula  [18] . 
The  formula  is. 


N+M 
=  Z 


x=l 


N+M 

Z 


y=l 


(2.16) 
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where 


Loss  formula  coefficients. 

Source  powers  (output  of  the  M  thermal 
and  N  hydroplants) 


2-5  System  Restrictions 


The  total  cost  of  operating  the  system  is  to  be 
minimized  subject  to  the  following  constraints: 

(i)  The  total  generation  of  the  system  should  be 
equal  to  the  system  load  demand  plus  the 
transmission  losses  in  the  system. 

That  is, 


M 

Z 

i=l 


Ti 


N 

Z  PHi 

j=i  1,3 


P  +  P 
D  L 


(2.17) 


This  equation  must  be  satisfied  at  all  instants 
during  the  optimization  interval. 

(ii)  The  amount  of  water  to  be  used  by  each  hydro¬ 
plant,  during  the  optimization  interval,  is 
predetermined  and  limited.  That  is, 

fT 

q.  dt  =  Q .  =  constant  (2.18) 

1  D 

Jo 

j  =  1 ,  2  ,  .  .  .  ,  N 

(iii)  Besides  the  above-mentioned  constraints,  the 

economic  schedule  must  conform  to  other  design 


•  T.  Gjs  03'  q 


16 


and  functional  requirements.  These  can  be 


classified  as, 


(a)  Design  limitations,  and 


(b)  Functional  limitations. 


The  design  limitations  appear  because  of  the 
reservoir  and  plant  characteristics.  The  size  of  the 
dam,  the  reservoir,  the  penstocks,  the  turbines,  etc., 
determine  these  limitations. 

The  limits  on  maximum  discharge  due  to  the  design 
considerations  can  be  written  as. 


i  qj 

Jmax  D 


(2.19) 


the  penstocks  and  the  turbines. 

The  forebay  elevation  is  restricted  to  a  definite 
range.  The  minimum  and  maximum  values  are  determined  by 
the  design  of  the  dam,  the  reservoir  and  the  location 
of  the  intakes.  That  is, 


(2.20) 


•  •  •  t 


N. 


The  forebay  elevation  is  directly  related  to  the 
storage,  hence  equation  (2.20)  can  also  be  interpreted  as. 
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s  . 

-’min  D 


<  s  . 
~  3 


—  s  • 

-’max  D 


j  =  1,  2, 


N 


(2.21) 


The  functional  limitations  arise  because  of  the 
fact  that  most  of  the  hydro-electric  plants  are  multi¬ 
purpose  in  nature;  hence  it  may  be  necessary  to  maintain 
certain  discharges  and  storage  levels  to  meet  obligations 
other  than  power  generation.  These  requirements  may  be 
seasonal  or  regular,  but  for  any  generation  schedule  to 
be  admissible,  these  conditions  must  be  satisfied.  Some 
of  these  restrictions  are  given  below. 

Maximum  forebay  elevation  due  to  flood  prospect, 


< 


max 


(flood) 


j 


=  1, 


N 


(2.22) 


Minimum  plant  discharge  and  spillage  to  meet 
the  irrigational  obligations. 


Jmin  (Irrigation) 
j  =  1,  2,  ....,  N 


(2.23) 


Minimum  plant  discharge  and  spillage  required  to  meet 
the  navigational  commitments. 


qj  +  °j  ~  qj  •  •  4--  \  (2.24) 

J  J  Jmin  (Navigation) 

The  storages  and  the  discharges  may  be  further 
restricted  by  project  limitations  other  than  those 


described  above. 


“xvrli  ed  VJ5  I  ;3  .dz  ill. 
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CHAPTER  III 

SCHEDULING  EQUATIONS  FOR  A  HYDRO-THERMAL 

SYSTEM  -  HYDRO- PLANTS  ON  SEPARATE  STREAMS 

3-1  Statement  of  The  Problem 

Consider  a  hydro-thermal  power  system  containing 
M  thermal  and  N  hydro-electric  power  plants,  interconnected 
electrically  by  a  network  of  transmission  lines.  The  hydro¬ 
electric  plants  are  located  on  separate  streams.  The  load 
demand  and  the  natural  inflows  are  known  functions  of  time. 
It  is  required  to  obtain  equations  for  computing  the  gene¬ 
ration  schedule  for  each  of  the  power  plants.  This  schedule 
must  be  such  that  the  operating  cost  of  the  system  over  the 
optimization  interval  is  minimized,  while  the  load  demand 
is  met.  The  transmission  losses  and  the  effect  of  varia¬ 
tion  of  head  upon  the  hydro-plant  characteristics  are  to 
be  taken  into  consideration.  Also,  the  system  restrictions 
explained  in  section  2-5  should  not  be  violated. 

3-2  Review  of  the  Current  Methods 

A  considerable  amount  of  work  has  been  done  to 
obtain  general  methods  and  equations  which  provide  the 
desired  economic  schedules.  A  large  number  of  these 
works  are  listed  in  references  [1,2]  of  this  thesis. 

Some  of  the  important  scheduling  equations  are  those 
derived  by  Ricard  [12],  Chandler,  Dandeno,  Glimn  and 
Kirchmayer  [13],  Cypser  [5],  Carey  [19],  Watchcorn  [4]  and 
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Glimn  and  Kirchmayer  [14].  It  has  been  shown  by  Glimn 
and  Kirchmayer  [14]  and  Arismunandar  [1]  that  these 
equations  are  essentially  equivalent  to  each  other. 

In  terms  of  the  symbols  adopted  in  this  thesis,  Glimn 
and  Kirchmayer  [14]  equations  are, 

th 

(i)  For  the  i  thermal  plant, 


9Ct  3P-r 

+  \  _J±  =  X 


(3.1) 


3  P 


Ti 


8PTi 


i  =  1 ,  2  ,  .  .  .  ,  M 


th 


(ii)  For  the  j  hydro-electric  plant, 


Yj0 


rt  3q  . 

dt 

Jo  3h  . 

3 

Aj 

3  q.  3Pt 

4d  l 

-  +  X  -  =  X  (3.2) 


3P 


Hj 


3P 


H  . 
: 


j  =  1 ,  2  ,  .  .  .  ,  N 


3-3  Mathematical  Formulation 


The  problem  stated  in  section  3-1  will  be  formulated 
in  the  following  paragraphs. 

The  cost  function  is  constructed  as  explained  in 
section  2-3.  The  time  integral  of  this  function  is  to 
be  minimized  over  the  optimization  interval,  i.e.. 


20 


I  = 


r T 
J  O 


PTM'  t}  dt  (3*3) 


is  to  be  minimized,  subject  to  the  following  constraints: 

(i)  The  total  generation  of  the  system  must  be  equal 

to  the  load  demand  plus  the  transmission  losses  in 
the  system,  i . e . , 


M 

Z 

i=l 


+  P, 


This  equation  can  be  rewritten  as, 


(3.4) 


M  N 

*  “  PD  +  PL  -  PTi  -  PHj  =  °  (3-5) 

1=1  j=l  ^ 

This  equation  must  be  satisfied  at  all  instants 
during  the  optimization  interval. 

(ii)  The  total  amount  of  water  to  be  used  by  each 

hydro-plant,  over  the  optimization  interval,  is 
predetermined  and  limited,  and  is  given  by, 


fT 
J  o 


dt 


constant 

j  =  1,  2,  ...,  N  (3.6) 
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(iii)  Various  design  and  functional  limitations 

(explained  in  detail  in  section  2-5)  when  com¬ 
bined  together,  result  in  putting  upper  and 
lower  limits  on  the  storages,  discharges  and 
the  power  outputs  of  the  hydro-electric  power 
plants.  These  can  be  expressed  as  follows: 
Discharge  limits: 


(3.7) 


represent  the  lower  and  upper  limits 


on 


Storage  limits: 

s  .  <  s  .  <  s’.  (3.8) 

-3-3-3 


j  =  1 ,  2  ,  .  .  .  ,N 

The  limits  on  storages  and  discharges  result  in 
limiting  the  power  output  of  the  hydro-electric  plants. 
These  can  be  expressed  as. 


(3.9) 


j  =  1 ,  2 ,  .  .  .  ,N 

Hence  I  in  (3.3)  is  to  be  minimized  subject  to 
the  constraints  given  by  equations  (3.5),  (3.6),  (3.7), 

(3.8)  and  (3.9).  The  two-sided  inequality  constraints 
given  by  (3. 7),  (3.8)  and  (3.9)  can  be  handled  by  transforming 


' 
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them  into  equality  constraints  [9] .  This  is  done  by 
introducing  a  new  variable  for  each  of  the  inequality 
constraints.  Such  manipulations  will  increase  the 
number  of  variables  of  the  system  considerably  (which 
ultimately  will  increase  the  number  of  Euler  equations 
for  the  system  by  the  same  number) .  For  this  reason 
these  constraints  are  not  considered  while  formulating 
the  problem. 

For  the  variational  formulation,  only  the  con¬ 
straint  given  by  equation  (3.5)  will  be  considered. 

The  minimization  of  I  subject  to  this  constraint  is  a 
Lagrange  problem  in  the  Calculus  of  Variations.  The 
problem  is  solved  by  replacing  the  cost  function  by 
an  augmented  function  H  [7,9],  defined  by. 


H  = 


A  <p 


(3.10) 


That  is. 


H  =  Cr 


P  +  P 
D  L 


M 

E 


N 


p  ,  -  y  p  . 
.  xi  ,VH] 
1=1  j=l 


(3.11) 


where  A  is  a  Lagrange  multiplier. 

Cm  is  given  by  equation  (2.15),  while  PT  and 

1  i-i 

P  .  are  defined  by  equations  (2.16)  and  (2.  7)  respec- 
tively.  Thus,  it  can  be  seen  that  H  is  a  function  of 
the  type , 


H  =  H  (P 


Tl 


T2 


,P 


tm' 


V 


V  V 


•  •  •  / 


SN'  X'  ^  (3*12) 


. 


. 
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Hence,  for  maximum  economy,  minimize 

fT 

J  ~  H(pT1'  PT2'  - 'PTM'  Sl'  S2  '  - 'SN' 

Sf'S  , . . . ,sN,X,  t)  dt  (3.13) 

The  extremal  arcs  are  obtained  by  solving  the 

Euler  equations  for  the  system.  Also,  for  any  extremal 

to  be  admissible,  the  constraints  given  by  equations 

(3.5),  (3.6),  (3.7),  (3.8)  and  (3.9)  must  be  satisfied. 

In  order  to  investigate  whether  the  functional  J  attained 

a  minimum,  it  would  be  desirable  to  check  that  the  other 

necessary  and  sufficient  conditions  are  satisfied  [1,  8]. 

However,  it  is  difficult  to  do  so  in  practice. 

Thus  the  scheduling  equations  for  the  power  plants 

are  given  by  the  Euler  equations  for  the  system. 

The  Euler  equations  for  the  system  are: 
th 

(i)  For  the  l  thermal  plant. 


3H 


3PTi 


=  0 


(3.14) 


i  =  1 ,  2  ,  .  .  .  ,  M 


th 

(ii)  For  the  j  hydro-electric  plant. 


3  H 


3  S  . 

D 


d  3  H 


=  0 


(3.15) 


dt  3s  . 
J 


j  =  1 ,  2  ,  .  .  .  ,  N 


„ 
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The  Euler  equation  for  the  hydro-electric  plants 
can  be  transformed  into  a  more  convenient  form,  as  shown 
in  the  following  section. 

3-4  Euler  Equation  in  the  Modified  Form 


The  coordinating  equations  for  the  hydro-electric 

plants  are  the  Euler  equations.  These  equations  are  put 

into  a  more  convenient  integro-dif ferential  form. 

fch. 

The  scheduling  equation  for  the  j  hydro-electric 
plant  is , 


3H 

3s  . 
3 

or 

3H 

3s  . 
3 

Let  us  define 


d 

3H 

dt 

3£  . 

3 

d 

3H 

dt 

3§  . 

3 

3H 

^  3H 

3s 

.  '  3S  . 

(3.16) 


(3.17) 


(3.18) 


only  when  _  is  not  equal  to  zero. 

3S  . 

3 

Note  that  in  this  thesis  it  is  assumed  that 


non-zero  at  all  times. 


is 


Hence 


3H  3H 


3  §  . 

3 


8  s  . 
3 


.  Z  (t) 


(3.19) 
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Now ,  from  equations  (3.17)  and  (3.19),  we  get 


3  H  d 

—  .  Z  (t)  =  — 

3  s .  J  dt 


3  H 


3  S 


or 


Z . (t)  dt 
D 


3H 


3  s  j 


3H 


3  S  . 

1  _ 


Integrating  both  sides,  we  get 


0 


Z . (t)  dt  =  log 
3  ^ 


3H 


3s 


+  K 


where,  is  an  arbitrary  constant  of  integration 
Let 


l°ge 

1  ' 

Y  . 

1  j0J 

3.22)  can  1 

then,  equation  (3.22)  can  be  rewritten  as. 


Z  .  (t)  dt  =  log 


0 


3  H 


Y 


jo 


3  s  . 


3J 


(3.20) 


3.21) 


(3.22) 


(3.23) 


(3.24) 
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or 


or 


1  3H 


jO 


3  H 


rT 

Z . (t)  dt 

Jo  3 

e 


rt 

Z  .  (t)  dt 
3 

J  0 

e 


(3.25) 


(3.26) 


where,  Z^ (t)  is  defined  by  equation  (3.18)  and  y j  q  is  an 
arbitrary  constant. 

This  equation  is  equivalent  to  the  Euler  equation 
and  will  be  called  the  "modified"  Euler  equation  in  the 
remainder  of  this  thesis. 


Derivation  of  the  Scheduling  Equations 

The  scheduling  equations  for  the  thermal  and  the 
hydro-electric  plants  of  the  system  are  given  by 
equations  (314)  and  (3.15)  respectively. 

th 

(i)  For  the  i  thermal  plant,  the  scheduling  equation  is, 

3  H 

-  =  0 


where  H  is  defined  by  equation  (3.11).  Hence,  we  get 


3P 


L 


3P 


Ti 


+  A 


A 


0 


(3.27) 
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or 


ac 


T 


3P, 


3PTi 


+  A 


=  A 


(3.28) 


ap 


Ti 


(ii)  For  the  hydro-electric  plants,  the  scheduling 
equations  are  obtained  by  using  the  "modified" 
Euler  equation  given  by  equation  (3.26) . 

t  h 

Hence,  the  scheduling  equation  for  the  j  hydro¬ 


electric  plant  is: 


aH 


as . 
3 


Z  .  (t)  dt 
0  3 


(3.29) 


jo  • 


where,  Z j  (t)  is  given  by  equation  (3.18). 


H 


is  defined  by  the  equation  (3.11),  while  PL 


and  P^j  are  given  by  equations  (2.16)  and  (2.7) 
respectively . 

Now  taking  the  partial  derivative  of  H  with 
respect  to  s^  and  substituting  into  equation  (3.26), 


we  aet 


3P, 


as . 

L  3 


ap 


Hi 


as . 

3  J 


=  Y 


jo 


Z . (t)  dt  (3.30) 


But , 


8PL  _  3PHj 


(3.31) 


as . 

3 


3PHj 


3Sj 


IO 


"bsilibo'i'i  orf*  pxtxcLi  Yd  lo  9XS  snoid 

. 


.tb  (i).S 


.9  sliriw  ,  (II.C)  noi&6£fp9  art  3  Yd  bani^sb 


■ 


. 
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and  putting, 


Z  .  (t)  dt 

0  3 


jO 


Yj(t) 


(3.32) 


in  equation  (3.30),  we  get. 


3P, 


L3PHj 


3PHi  .  ffHl 


3  S  . 

3 


3s 


3  J 


=  Y j  (t) 


(3.33) 


or 


3P 


Hi 


3s  . 
3 


3P, 


-  1 


3P 


Hj 


=  Y j  (t) 


(3.34) 


or 


3P, 


aPHj 


-  1 


=  Y  •  (t)  .  Tp - 

^  !!si 


3s  . 
3 


(3.35) 


Now,  consider 


3s  .  1 

— =  ap 

3P„_.  Hj 


Hj 


3s 


(3.36) 


3P 

This  equation  is  true  if  the  derivatives  Hj  and 

3S  . 

3 s  ■  exist  and  are  non-zero.  A  proof  is  given  ^  by 
8PHj 

Glimn  and  Kirchmayer  [14] . 


i  -■  v  .r.  ..  I  o,  i  '• 
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Hence,  equation  (3.35)  can  be  rewritten  as. 


3P, 


L8PHj 


-  1 


ds . 

=  Y  •  (t)  .  - 

3  3PHj 


(3.37) 


Rearranging  the  terms ,  we  get 


Yj(t)  . 


3s  . 

.  -J- 

9P«j  J 


3P, 


+  X 


3PHj 


=  A 


(3.38) 


Y . (t)  is  given  by  equation  (3.32),  while  Z . (t) 

J  i  ^ 

is  given  by  equation  (3.18)  as, 

3H 
3S  . 

Z  .  (t)  =  - 3- 

J  3H 

3  s  . 

3 

3  H 

In  order  to  evaluate  Z  .  (t) ,  the  values  of  - —  and 

3  3s . 

9  H  ^ 

-r-r—  are  to  be  obtained. 

9  S  . 

3 

Taking  the  partial  derivative  of  H  with  respect 


to  s j ,  we  get 


3H 


3s  . 
3 


=  X 


3P, 


3  S  . 

[_  3 


8PHj 

3  S  . 

3 


(3.39) 


but. 


3P, 


3s 


3P, 


3P 


Hj 


!!si 

3S  . 

3 


(3.40) 


IW  2Ti  ■  r  --1-1  : : 
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hence 


9H 


9s  . 

3 


(3.41) 


or 


9H 


9s 


=  X 


!!hi 

9S  . 

3 


9P, 


LaPHj 


-  1 


(3.42) 


Now,  taking  the  partial  derivative  of  H  with 
respect  to  § ^ ,  we  obtain 


9H 


9  S  . 
3 


=  X 


9P, 


9  S  . 

L  3 


3PHi 

9  s  . 

3  _ 


(3.43) 


9P 


Substituting  for  L  from  equation  (3.31)  in 

9s  . 

3 

equation  (3.43)  and  rearranging  the  terms,  we  get 


9H 


9s 


=  X 


9PU  . 
--»! 

9s  . 

3 


9P, 


_8PHj 


-  1 


(3.44) 


Now,  from  equations  (3.42),  (3.44)  and  (3.18) 

we  obtain 


x  3P«i 

PPl  -  il 

9  S  . 

3 

K  J 

x  3PHj 

PPl  -  il 

9s  . 

3 

8PHj 

(3.45) 


. 


' 

jag  sw  v8 mi9i  drf*  bns  (£&.£)  fio.< 
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or 


Zj  (t) 


(3.46) 


Hence  from  equations  (3.38),  (3.32)  and  (3.46), 

th 

we  get  the  scheduling  equation  for  the  j  hydro-electric 
plant  as, 


•  t 

0 


a  s . 

L  3 


dt 


3s  . 

3 


8PHj 


3P. 


+  X- 


spHj 


X 


(3.47) 


3-6  Comparison  with  the  Scheduling  Equations  Derived 
by  Glimn  and  Kirchmayer  tl4T~ 

th 

Their  scheduling  equation  for  the  l  thermal 

plant  is  given  by  equation  (3.1).  This  is  identical  to 

that  derived  in  this  thesis  (given  by  equation  (3.28)). 

th 

For  the  j  hydro-electric  plant,  their  scheduling 
equation  is  given  by  equation  (3.2). 


. 


Also,  they  have  shown  in  their  paper  [14]  that 


3P. 


M 


3h 


3P 


9qj 


(3.48) 


Thus,  equation  (3.2)  can  be  rewritten  as. 


Yj0  e 


0 


9PHi  9PHj 


3h  . 
3 


3q. 


dt 

Aj 


3q  . 
~3 


9PHj 


3P, 


+  X 


8PHj 


=  A 


(3.49) 


Now,  comparing  the  equations  (3.47)  and  (3.49),  we  find 
that  the  equations  are  identical  only  if  the  following 
relations  are  true: 


3s  .  3q . 
-  3  =  3 


3P 


Hj 


3P 


Hj 


(3.50) 


and , 


(3.51) 


The  discharge  of  a  plant  can  be  represented  as  the 
difference  between  the  natural  inflow  and  the  rate  of 
change  of  storage  in  the  plant's  reservoir,  i.e.. 


(3.52) 


also,  we  know  that, 


. 


' 

,  P€ 

.  rid.i  biz 6  wo  In  :  iu  3  1  i  ‘ 


33 


hence,  it  can  be  concluded  that  both  q.  and  (F  -  s.) 

3  3 


are  functions  of  PTT .  and  h . 

H3  3 


F .  is  an  uncontrollable  variable 
3 


and  is  also  independent  of  j .  By  taking  the  partial  deri¬ 
vatives  of  both  sides  of  equation  (3.52),  we  get 


8q 


2 


as 


(3.54) 


ap 


Hj 


ap 


Hj 


Hence,  equation  (3.50)  is  true. 

Furthermore,  if  the  tail  race  elevation  and  the 

losses  in  the  conduits  are  taken  to  be  constant,  and  if 

th. 

the  surface  area  of  the  reservoir  of  the  j  hydro¬ 
electric  plant,  A_.  ,  is  a  constant,  then  the  net  head 
at  the  plant  can  be  written  as. 


h .  =  sj  +  K,  . 

3  aT"  ^ 

3 


(3.55) 


where  K,  .  is  a  constant. 
13 


Again , 


Therefore 


PHj  PHj  (qj'  hj'  t) 


=  PHj  [V  (4r  +  Klj>'  (3‘56) 


aPH-i  8PHi 

— £2  =  A  — 22 


ah. 

3 


a  s  j 


(3.57) 
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or 


9P 


A  . 
3 


HI  = 


9P 


ah 


HI 


as  . 

3 


(3 


Therefore,  equation  (3.51)  is  true. 

Thus,  the  scheduling  equation  (3.47),  derived  in  this 

th 

thesis  for  the  j  hydro-plant,  is  equivalent  to  that 
of  Glimn  and  Kirchmayer  [14] . 


58) 
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CHAPTER  IV 

SCHEDULING  EQUATIONS  FOR  A  HYDRO-THERMAL 

SYSTEM  -  SOME  OF  THE  HYDRO-PLANTS 
LOCATED  ON  COMMON  STREAMS 

4-1  Statement  of  the  Problem 

The  power  system  under  study  consists  of  M  thermal 
and  N  hydro-electric  power  plants,  interconnected  electri¬ 
cally  by  a  net-work  of  transmission  lines.  For  simplicity, 
it  is  assumed  that  all  except  two  hydro-plants  are  located 
on  separate  streams.  Let  the  k  and  k+l°  be  the  two 
hydro-plants  located  on  the  same  stream.  We  assume  that 
the  load  demand  and  the  natural  inflows  are  known  functions 
of  time.  The  time  taken  by  the  water  to  flow  from  the 
upstream  plant  to  the  downstream  plant  is  t.  It  is 
desired  to  obtain  the  scheduling  equations,  such  that 
the  most  economical  generation  schedule  can  be  obtained. 

4-2  Review  of  Current  Methods 

In  1941,  Burr  [20]  made  an  attempt  to  develop 
a  general  method  for  determining  the  most  economical 
loading  of  the  common-flow  hydro-electric  plants. 

Using  the  Steepest  Descent  method,  Cypser  [5]  solved 
the  Euler  equations  for  a  system  in  which  the  two  hydro¬ 
electric  plants  are  located  on  the  same  stream,  Menon 
[15]  solved  the  common-flow  problem  by  constructing  a 


' 


1 
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series  of  minimizing  sequences.  Kirchmayer  [18] 
extended  the  incremental  cost  technique  to  the  hydro¬ 
plants  located  on  common  water-shed.  The  hydro-plants 
were  treated  as  constant  head  plants.  Drake,  Kirchmayer, 
Mayall  and  Wood  [21]  used  that  method  to  coordinate  the 
expenditure  of  the  steam  and  hydro-electric  resources  of 
the  South  California  Edision  Company.  The  effect  of 
variation  of  head  upon  the  hydro-plant  characteristics 
was  not  taken  into  consideration. 


4-3  Mathematical  Formulation 


The  problem  stated  in  section  4-1  is  similar  in 
nature  to  that  stated  in  section  3-1.  The  only  difference 
is  that  the  k  and  k+lb  hydro-plants  are  located  on 
the  same  stream.  The  following  analysis  is  quite  general 
in  nature  and  can  be  extended  to  a  system  in  which  more 
than  two  hydro-plants  may  be  located  on  the  same  stream. 

The  cost  function  for  the  system  is  constructed 
as  explained  in  section  2-3.  The  time  integral  of  this 
function  is  to  be  minimized  over  the  optimization 
interval.  That  is,  we  minimize 


I 


PTM'  dt  (4,1) 


subject  to  the  constraints  listed  in  section  2-5. 


[ 8X1  TaYfifwto*** 

...0-  fc*  o 

.  IC  s  CEJ  alni  ns:  s; 
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The  only  constraint  which  will  be  considered  here  is 
that  the  total  generation  of  the  system  must  be  equal 
to  the  load  demand  plus  the  transmission  losses  in  the 
system,  that  is, 

M  N 

i=l  PTi+j^  PHj  =  PD  +  PL  (4'2) 


or 

*  =  PD  +  PL  -  E  0  (4‘3) 

However,  for  any  optimal  solution  to  be  admissible 
all  other  constraints  described  in  section  2-5  must  be 
taken  into  account. 

Again,  minimization  of  I  in  (4.1)  subject  to  the 
constraint  given  by  equation  (4.3),  is  a  Lagrange  problem 
in  the  Calculus  of  Variations.  It  is  handled  by  replacing 
CT  by  an  augmented  function  H  [7,9],  defined  by, 

H  =  CT  +  X  <f> 

M  N 

=  CT  +  MPD+PL  -  (4‘4) 

From  equations  (2.16),  (2.7)  and  (2.11),  it  is 

seen  that  H  is  a  function  of  the  type, 

H  =  H  (PT1'  PT2 '  ‘ '  PTM'  Sl'  S2 '  SN' 

Sl'  S2  ’  SN'X't)  (4,5) 


.ir/pa  d  '  jm  msdwa  adi  *  no:rJs-  p  J  1 


dsi/pa  y<*  navip  inir  ii3/ioo 
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Now,  for  maximum  eqonomy,  we  minimize 


(4.6) 


The  optimal  solution  is  obtained  by  solving  the 
Euler  equations  for  the  system.  The  Euler  equations  for 
the  system  are, 


3H 

3P 


Ti 


0 


i  «  1,  2, 


N 


(4.7) 


and 


3H  d 

3  s  .  dt 
3 


0 


j  =  1,  2, 


(4.8) 


N 


Equation  (4.8)  can  be  transformed  into  a  more 
convenient  form  as  given  by  equation  (3.26)  and  will  be 
used  for  obtaining  the  required  scheduling  equations. 


4-4  Derivation  of  the  Scheduling  Equations 

The  scheduling  equations  for  thermal  and  the 
hydro-electric  plants  are  obtained  by  using  equations 
(4.7)  and  (3.26)  respectively. 


■ 
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(a)  Taking  the  partial  derivative  of  H  (defined  by 
equation  (4.4))  with  respect  to  PT^  and  substi¬ 
tuting  into  equation  (4.7),  we  get, 


dCr 


3P, 


3  P  . 

Tl 


+  A 


3  P_  • 
Tl 


-  A  =  0 


or 


3Cr 


3P. 


3PTi 


+  A 


3PTi 


=  A 


(4.9) 


(4.10) 


This  equation  is  identical  to  that  for  the 
th 

i  thermal  plant  of  the  system  discussed  in  Chapter- 

3. 

(b)  The  scheduling  equations  for  the  hydro-plants 

are  obtained  by  using  the "modified"  Euler  equation 
(3.26)  . 

(i)  The  hydro-electric  plants  situated  on 
separate  streams: 

The  coordinating  equations  can  be  derived  as 

in  section  3-5.  The  coordinating  equation  for  the 
th 

j  hydro-plant  is  (refer  to  equation  (3 . 47) )  . 


Y 


j0 


3P 


3  S  . 

L  3 


/ 


dt 


+  A 


3P 


Hj 


A 


(4.11) 


j  =  lf  2,  k-l,k+2,...,  N 

j  ^  k,  k+1 
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s  t 

(ii)  The  downstream,  k+1  ,  hydro-electric  plant: 

The  coordinating  equation  is  obtained  by 
taking  the  partial  derivatives  of  H  with 
respect  to  s^.+^  and  and  substituting 

them  in  the  ifiodified  Euler  equation  (3.26). 
Hence,  the  scheduling  equation  for  the 
hydro-plant  is. 


3P, 


3P 


Hk+1 


3S 


k+1 


3  S 


k+1 


0 


zk+i<t)  at 


=  Y 


(k+1) 0  * 


(4.12) 


where  Z^.  (t)  is  given  by  equation  (3.18) 
Also , 


3P, 


3P. 


3  P 


Hk+1 


3sk+l  3PHk+l  3  sk+l 


(4.13) 


because  appears  in  the  expression  for  P^+2  only 

(refer  to  equation  (2.11)). 


9  P 

Substituting  for  L 


9  s 


k+1 


in  equation  (4.12)  we  get 


3  P, 


3  PHk+l  3PHk+l 


3PHk+l  3  Sk+1  3sk+l 


zk+i(t)  dt 


=  Y 


(k+1) 0 


.  e 


(4.14) 


,-j  • :  1  4  ?.  O  a  >v  '  ; 


' 


• 
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or 


3P 


Hk+1 


3  S, 


3P, 


-  1 


3  P 


Hk+1 


,0 

Y  (k+1) 0  6 


Zk+l(t)  dt 


(4.15) 


also,  it  can  be  shown  [14]  that 


3  P 


Hk+1 


=  3§ 


3  s 


k+1 


k+1 


3  P 


Hk+1 


(4.16) 


Hence,  equation  (4.15)  can  be  rewritten  as, 


(k+1) 0 


.  e 


zk+i(t)  dt 


3S 


k+1 


3P 


Hk+1 


3P, 


+  A 


=  A  (4 


3  P 


Hk+1 


In  order  to  find  the  expression  for  Z^+^(t) ,  the 
partial  derivatives  of  H  with  respect  to  s^+^  and  s^+^ 
are  to  be  evaluated. 


but , 


3  H 


=  A 


3  S 


k+1 


3  P, 


3  P 


Hk+1 


3  S 


k+1 


3  S 


k  +  1 


3P, 


3P, 


3  P 


Hk+1 


3sk+l  3PHk+l  3sk+l 


(4.18) 


(4.19) 
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Hence,  substitutina  (4.19)  into  (4.17)  we  aet 


3H 


3P 


=  X 


Hk+1 


3s 


k+1 


3  S 


k+1 


3P, 


3s 


k+1 


Similarly,  we  obtain 


-  1 


(4.20) 


3H 


3  P 


=  X 


Hk+1 


3  s 


k+1 


9  S 


k+1 


3P. 


3  S 


k+1 


Substituting  for 


3  H 


we  aet 


Zk+1 (t) 


3s 


and 


-  1 


3H 


(4.21) 


k+1 


3S 


in  equation  (3.18), 


k+1 


3P 


Hk+1 


3  S 


k  +1 


3P, 


-  1 


3  P 


Hk+1 


3P 


Hk+1 


3  S 


k  +1 


3  P, 


-  1 


3  P 


Hk  +1 


(4.22) 


or 


zk+i(t) 


3P 


Hk+1 


3  S 


k+1 


3P 


Hk+1 


3  S 


k+  I 


(4.23) 


Hence,  from  equations  (4.17)  and  (4.23),  the 


scheduling  equation  for  the  k+1  hydro-plant  is 

3P, 


Y 


(k+1) 0 


.  e 


•  t 

9PHk+l 

J  0 

—  sk+l 

/ 


st 

(Hk+l 


9  s 


k+1 


dt 


3S 


k+1 


3P 


Hk+1 


3P, 


+  X- 


=  X 


3  P 


Hk+1 


(4.24) 


u 


, 
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This  equation  is  identical  to  the  scheduling  equation 
for  the  hydro-plants  located  on  separate  streams  (equation 
(4.11)  )  . 


'tl'l 

(iii)  The  upstream,  k  (  hydro-electric  plant: 

The  water  released  from  the  upstream  plant 
at  an  instant  (t-x)  reaches  the  downstream  plant  at  an 
instant  t.  This  discharge  flows  into  the  reservoir  of 
the  downstream  plant  and  thus  affects  the  operation  of 
that  plant.  This  has  been  discussed  in  detail  in 
section  2-2.  There  it  has  been  shown  that  the  variable 

.  P  p 

s^  appears  in  the  expressions  for  both  Fk  and  Hk+1 
(refer  to  equations  (2.8)  and  (2.11)),  while  the  variable 

s^.  appears  only  in  P  ^ .  This  must  be  kept  in  mind  when 

taking  the  partial  derivatives  of  H  with  respect  to  s^ 

and  s,  . 
k 


.  th 

The  scheduling  equation  for  the  k  plant  is 
obtained  by  taking  the  partial  derivatives  of  H  with 
respect  to  s^.  and  s^. ,  and  substituting  them  into  the 
''modified'5  Euler  equation  (equation  (3.26)). 

Thus ,  we  get 


9P, 


9  P 


Hk 


9  P 


Hk+1 


9  S, 


9  s, 


9S, 


=  Y 


kO 


0 


Z.  (t)  dt 
k 


(4.25) 


where,  (t)  is  given  by  equation  (3.18). 
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Also,  we  can  write 


3P, 


3P, 


3PHk  .  3PL 
— -  - - 


3P 


Hk+1 


3  S,  3P„,  3  S. 

k  Hk  k 


3PHk+l  9sk 


because  P^  is  a  function  of  both  P  ^  and  while 

both  P  ^  and  P^+^  are  functions  of  s^. 

Hence,  equation  (4.25)  can  be  rewritten  as. 


3P, 


BPHk  !!l 


9PHk  9  sk 


3P 


Hk+1 


3P 


Hk+1 


3P 


Hk 


3  P 


Hk+1 


3  s, 


3  s 


3  s. 


kO 


k 

ft 

lo  Zk(t)  dt 


.  e 


Now,  putting 
t 


Y 


0 


Z,  (t)  dt 
k 


kO 


=  ^k(t) 


and  rearranging  the  terms,  we  get 


3  P 


Hk 


3  s. 


3P, 


-  1 


3  P 


Hk 


3  P 


+  X 


Hk+1 


3  s. 


3P, 


-  1 


3  P 


Hk+1 


=  Yk(t) 


or 


3P 


Hk 


3  Sn 


3P, 


—  PHk 


-  1 


3P 


=  Yk  (t)  -  X 


Hk+1 


3  s. 


3P, 


-1 


3  P 


Hk+1 


4.26) 


4.27) 


4.28) 


4.29) 


4.30) 


- 


' 


. 
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We  can  write  [14]  , 


3P 


Hk 


3  s. 


3S, 


3P 


Hk 


Hence,  equation  (4.30)  can  be  written  as, 


3P, 


3PH, 


-  1 


3  P 


Yk(t)  ~  X 


Hk+1 


3S, 


3P, 


-1 


3P 


Hk+1 


3  S, 


3P 


Hk 


or 


3P 


Yk(t)  "  X 


Hk+1 


3S, 


3P, 


3P 


-1 


Hk+1 


3S, 


3P 


Hk 


■h  X 


*  X 


3  P 


HL 


Rearranging  the  terms  of  the  scheduling  equation 
s  t 

for  the  k+1  hydro-plant  (equation  (4.24)),  we  get 


3P, 


-  1 


3PH 


k+1 


3s 


=  Yk+i(t)  • 


k+1 


3P 


Hk+1 


where 


t 

0 


Yk+l(t)  “  Y  (k+1 ) 0  * e 


3  P 


Hk+1 


/ 


3s 


k+1 


3  P 


3  S 


Hk+1 


k+1 


dt 


(4.31) 


(4.32) 


(4.33) 


(4.34) 


(4.35) 


A  ~ 


-4 
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Now  substituting  for  A 


(4 . 33) ,  we  get , 


3P, 


3  P 


-  1 


Hk+1 


in  equation 


3  S 


Yk(t)  ‘  Yk+i(t)  • 


k+1 


3  P 


Hk+1 


8PHk+l  35k 


3s,  -t 
k 


3P 


Hk 


3P, 


+  A 


=  A 


3P 


Hk 


(4.36) 


y^(t)  and  y^+^(t)  are  given  by  equation  (4.28)  and 
(4.3  5)  respectively;  and  Z^(t)  is  given  by  equation  (3.18) 


In  order  to  evaluate  Z^(t)  we  must  find  the  partial 
derivatives  of  H  with  respect  to  s^  and  s^. 

The  variable  s^  appears  only  in  P  ^  /  hence,  by 
taking  the  partial  derivative  of  H  with  respect  to 
s^.,  we  get, 


3  H 


3  s, 


=  A 


3  P, 


-3sk 


3P 


Hk 


3  S, 


(4.37) 


Also 


3  P  3Pt 

i_j  1j 


therefore 


3P 


Hk 


3sk  3PHk  3sk 


(4.38) 


3  H 


3  s. 


=  A 


3P, 


3  P 


Hk 


3  PHk  3sk 


3  P 


Hk 


3  Sn 


(4.39) 
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or. 


8  H 


3  S, 


3P 


Hk 


3  S, 


3  P, 


-  1 


3  P 


Hk 


(4.40) 


Since  both  P^  and  are  functions  of  ,  we 


obtain , 


3  H 


3  S, 


=  X 


9PL  _  8PHk  _  9PHk+l 
3sk  3sk  3sk 


(4.41) 


and  using  equation  (4.26) ,  we  can  write. 


3  H 


3  s, 


=  X 


3  P, 


3PHk  3PL 
— - -  +  - 


9PHk+l  3PHk 


3  P 


Hk+1 


3PHk  9Sk 


9  PHk+l  3*k 


3  s, 


3  S, 


(4 


rearranging  the  terms,  we  get 


3  H 


3  s, 


3  P 


=  X 


Hk 


3s, 


ph  .p 

9PHk+l 

1PL  _  , 

_?PHk 

r  A 

3  S, 

k 

_9PHk+l 

(4.43) 


3  H  3  H 

Now,  substitutina  for  - — —  and  — r—  in  equation  (3.18), 

3sk  3§k 


we  get 


3P 


Hk 


3  S 


k 


zk(t)= 


3  P, 


-  1 


3  P 


Hk 


3  P 


Hk 


3  S, 


3  P, 


-  1 


3P 


Hk 


3  P 


+  X 


Hk+1 


3  S, 


(4.44) 


3  P, 


-  1 


3  P 


Hk+1 


42) 


(I*  M 


<£*.M 
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or 


Zk(t)  = 


9P 


Hk 


3P. 


-  1 


3sk 

313 

—  Hk 

SPHk 

PpL  'J 

3PHk+l 

PPl  -  il 

3sk 

_3PHk  _ 

i 

as. 

k 

_9PHk+l 

(4.45) 


th 

Thus,  the  scheduling  equation  for  the  k  plant  is 
given  by  equation  (4.36)  ,  where  Y^(t)  and  Yjc+^(t)  are  given 
by  equations  (4.28)  and  (4.35)  respectively.  Furthermore, 
Z^(t)  is  defined  by  equation  (4.45). 


4-5  Discussion  of  the  Scheduling  Equations  for  the 
Common-flow  Plants 


If  we  substitute 


Yk(t) 


Yk+1 (t) 


3sk+l  3PHk+l  =  Yk 
3PHk+l  3sk 


(4.46) 


in  equation  (4.36),  we  get 


Y  i 


3  S, 


3Pul 

L.  Hk  j 


SP, 


+  X 


=  X 


9P 


Hk 


(4.47) 


Comparing  this  equation  with  equation  (4.24)  for 
s  t 

the  k+1  hydro-plant,  it  is  noticed  that  both  the  equa¬ 
tions  are  similar  in  nature.  The  difference  lies  in  the 


way  and  yk+1 


vary  with  time. 
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th 

Y j  for  the  j  (j^k)  hydro-plant  is  interpreted 
as  the  water  conversion  coefficient  for  the  plant  [14] . 
Its  value  determines  the  amount  of  water  which  is  used 
over  the  optimization  period.  The  influence  of  these 
coefficients  on  the  scheduling  equation  is  discussed  in 
section  5-3.  The  water  conversion  coefficient  for  the  up- 

s  t 

stream  plant  is  modified  due  to  the  presence  of  the  k+1 
(downstream)  plant  on  the  same  stream.  Thus  if  the 
value  f°r  I  (k+i)  o  -*-s  Ranged,  the  value  for  is  also 
affected. 

The  time  lag  between  the  two  reservoirs,  x,  is 
involved  in  the  expression  for  due  to  s^(t-x). 

In  the  scheduling  equation,  x  affects  only  the  value  of 
3P, 


Hk+1  ,  and  this  appears  only  in  the  coordinating 


9  S, 


th 


equations  for  the  k  hydro-electric  plant. 
Now,  from  equation  (2.9),  we  have 


PHk+l(t)  PHk+l  (sk+l(t)/  (sk(t  x)  +  sk+l(t))) 


and  as  the  variables  of  the  system  can  be  varied  indepen¬ 
dently  , 


3PHk+l (t)  _  3PHk+l (t) 

a^k+l(t)  9sk(t-x) 


£  ?;  ex  rfoxxiw  z  5dj  w  do  ini 

^X+d  9riJ  do. 


L+3tHq 


-  *q  fti  •  >xxrv  ed  rxso  -Y  add  do  a  alt'.  ■/  5 
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Also,  s^(t-x)  can  be  written  as  some  non-linear 
function  of  (t)  and  t,  i.e., 

sk (t-x )  ®  f  (sk (t) ,x ) 

Then 

3PHk+l (t)  _  3 PHk+l 
3sk(t)  3sk(t-r) 

Using  equation  (4.48)  ,  equation  (4.50)  can  be  written  as, 

3PHk+l  3PHk+l(t>  .  d®k(t-T) 

- -  =  — - -  — - -  (4.1 

3  ^  8Sk+i(t)  dsk(t) 


3Sk  (t-T ) 


3sk(t) 


(4.49) 


(4.50) 


Now,  if  t  is  constant,  sk(t-x)  can  be  expanded  using 
a  Taylor's  series,  as 


sk(t-x)  =  sk(t)  -  x  sk(t)  +  x_^  s’  (t)  - 


(4.52) 


Differentiating  both  the  sides  with  respect  to 


sk (t) ,  we  get 


dsk(t-x)  ds*k  (t)  x  d*S k  ( t ) 

- ; -  =  1  -  X  - : -  +  -  - 

dsk (t) 


dsk(t)  2!  dsk(t) 


Hence,  a  suitable  numerical  method  for  evaluating 

etc.  must  be  evolved  in  order  to  obtain  the  value 
3s,  (t-x) 


(4.53) 


dsk(t) 
dsk (t) 


of 


th 


3sk (t) 


.  Then,  the  scheduling  equation  for  the 


k  hydro-plant  can  be  solved  numerically. 


eul&v,  ©rid  rciRddo  o:t 
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CHAPTER  V 
APPLICATION 

5-1  The  Problem  Considered 

The  system  considered  here  consists  of  one  thermal 
and  two  hydro-electric  plants.  It  is  required  to  operate 
this  system  to  meet  a  given  load  demand  in  the  most 
economical  manner.  The  characteristics  of  the  hydro-plants 
and  the  cost  function  for  the  system  are  given  in  Appendix 

1.  The  system  is  shown  schematically  in  Fig.  5.1  The 
optimization  interval  is  24  hours.  The  two  hydro-plants 
are  located  on  the  same >stream.  The  effect  of  variation 
of  head  upon  the  hydro-plant  characteristics  will  be  taken 
into  consideration. 

The  problem  is  solved  under  the  following  assumptions 

1.  The  time  taken  by  the  water  to  flow  from  the  upstream 
plant  to  the  downstream  plant  is  assumed  to  be 

zero . 

2.  During  the  optimization  period,  the  natural  inflows 
to  the  reservoirs  of  the  two  plants  are  assumed  to  be 
zero . 

3.  The  tail  race  elevations  and  the  head  losses  in 
conduits  are  assumed  to  be  constant. 

4.  The  transmission  losses  in  the  system  are  negligible. 


i 

senoxiqr  ;i  Coi  orf^  iQhnu  b9vIoa  si  meJ 
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System  Studied  in  Chapter  5 


Transmission  Lines 


FIGURE  5-1 


53 


5-2  Method  of  Solution 


The  above-mentioned  problem  is  solved  by  using 
the  scheduling  equations  developed  in  section  4-4. 

The  equations  are  simplified  in  accordance  with  the 
assumptions  listed  in  the  previous  section. 

Neglecting  the  terms  containing  P  in  equation 
(4.10),  we  get  the  scheduling  equation  for  the  thermal 
plant  as, 


3C 


T 


3P 


T 


=  X 


(5.1) 


Similarly  the  scheduling  equation  for  the  down¬ 
stream  plant  (plant  no.  2)  becomes  (equation  4.24)): 


-  y(2)0 


•  t 

9PH2  ,  3PH2 

*  • 

Jo 

3s2  3s2 

at 


X 


(5.2) 


and  from  equation  (4.36),  we  obtain  the  scheduling  equation 
for  the  upstream  plant  (plant  no.  1)  as. 


Y 1  ( t )  -  Y2  (t) 


3  s. 


3P 


H2 


3pH2  9s1 


3s. 


3P 


H2 


=  X 


(5.3) 


X 
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where 


Yl(t)  =Y(l)0-e 


and 


Y2(t)  Y  ( 2 )  0  ’ e 


ap 


HI 


8  s  • 


3PH1  +  3PH2 


as. 


as. 


dt 


ap 


H2 


as. 


ap 


H2 


as. 


dt 


(5.4) 


(5.5) 


Equation  (5.3)  can  be  further  simplified  by  taking 
into  account  the  fact  that  t  =  0.  When  t  =  0 ,  equation 
(2.9)  becomes, 

PH2  =  PH2  ^ s2 '  +  ®2*'  ^  (5.6) 

By  taking  the  partial  derivatives  of  both  the  sides 
it  is  found  that 


3PH2  3PH2 


as. 


as. 


Also,  we  know  that  [14], 


as. 


ap 


H2 


ap 


H2 


as. 


(5.7) 


(5.8) 


■ 
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therefore 


3  s. 


3  P 


H2 


=  1 


(5.9) 


3pH2  Zs± 


Hence, equation  (5.3)  can  be  rewritten  as, 

3  s, 


Yf  ( t )  "  Y2(t) 


( 


3  P 


-)  =  A 


(5.10) 


H2 


Equations  (5.2)  and  (5.10)  can  be  rewritten  as, 


3P 


*  Y2(t)  = 


H2 


(5.11) 


3s. 


and 


1 

A 


Y 1 ( t )  ~  Y2(fc) 


(5.12) 


respectively.  (t)  and  Y2  (t)  are  given  by  equations 
(5.4)  and  (5.5)  respectively. 

Now,  from  the  given  system  data  (Appendix  1) 

CT  =  0.012  PT2  +  4.0Pt  +  1.0  (5.13) 

Hence,  evaluating  the  partial  derivative  of  CT  with 
respect  to  PT  and  substituting  into  equation  (5.1),  we  get 

0.024  PT  +  4.0  =  A  (5.14) 


or 


A  -  4.0 
0.024 


(5.15) 


* 
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Equation  (5.15)  represents  the  scheduling  equation 
for  the  thermal  plant. 

Also,  from  the  system  data  (Appendix  1)  for  the  down¬ 
stream  plant  (plant  No.  2),  we  have. 


PH2  =  0.076  q2  (h2  -  5  -  1.5  q2) 


(5.16) 


Now,  since  the  natural  inflows  are  zero  equation 
(2.5)  becomes 


q2  =  -Sl  -  S2 


(5.17) 


Also 


2  ^ 


(5.18) 


Substituting  for  q9  and  h2  inequation  (5.11),  we  get 

s~ 


PH2  =  °*076  (~si_s2 ) 


A  "  5  -1. 5  (  s1  s2) 

2 


(5.19) 


Hence , 


3PH2  =  3PH2 
3s^  3s2 


=  -0.076 


(5.20) 


-  5  -  3.0  (-s1-s2) 


=  -0.076  (h2  -  5  -  3.0  q2) 


(5.21) 


(,s-  i-)  e.x 
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and 


3P 


H2 


3  s. 


0.076  -=r=- 

A2 


Using  equations  (5.21  and  (5.22),  equation  (5.2) 
for  the  downstream  plant  can  be  written  as. 


Y2(t) 


=  0.076  (h,  -  5  -  3.0  q,) 


and,  from  equation  (5.5) 


•t  0.076  q2 


0 


-A2. 0. 076 (h2~5. 0-3.0  q2) 


dt 


Y2 (t)  Y (2) 0  ,e 


Y (2) O.e 


-  ^2 

0  A2 . (h2-5 .0-3.0  q2) 


dt 


Thus  , 


Z2(t)  = 


A2  (h2~5 .0-3.0  q2 ) 


Now,  by  solving  equation  (5.23)  for  q2 ,  we  get. 


q2 


1 

3 


Y2 

0.076.X 


-  h2  +  5.0 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


(5.26) 


(£.£)  nox^Bupa’ , (SS.2 
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From  the  system  data  (Appendix  1)  the  expression  for 
the  output  of  the  upstream  plant  (plant  No.  1)  is, 

PH1  =  0. 06486. gx  (hx-20  +  38.107  -  2.863  q^2)  (5 

Since,  the  natural  inflows  are  assumed  to  be  zero, 


q 


1 


(5.28) 


Also 


h 


1 


(5.29) 


Substituting  for  and  h^  in  equation  (5.27)  and 
taking  the  partial  derivatives  with  respect  to  s^  and  s^, 
we  get 


and 


3P 


HI 


3s. 


0.06486  (-§  ) 
A, 


0.06486  q^ 


(5.30) 


3  P 


3  S 


—  =  -0.06486  (h^-20+76. 214  (-s.^-8. 589  (-s^2) 


=  -0.06486  (h1-20+76. 214  q±  -  8.589q1  ) 
Substituting  for  3PH1  in  equation  (5.12),  we  get 


3  s. 


.27) 


(5.31) 


'  3  ■{'  •  : 


’ 
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1 

X 


+  0.06486  (h^-2 0+76. 214  q 


Rearranging  the  terms  and  solving  for 


76.214  +  /  ( 76 . 214 )  2  -  34.356  C1 

17.178 


where 


Yl“  Y2 


X*  0 . 06486 


hl  +  20 


Equation  (5.6)  can  be  written  as 

t 


( t)  dt 


Y1  ^  Y  (1)  0  e  0 


where 


ap 


HI 


as. 


z1(t) 


ap 


HI 


ap 


H2 


as. 


as. 


Substituting  for  3Phi  3Phi  and  3Ph2 


as1  '  asx 


a T 


(5.30),  (5.31)  and  (5.21)  in  equation  (5.36) 


-  0.06486  .  q^ 


z1(t)  = 


0.06486  (h1-20  +  76. 214  q 


+  0.076  (h2-5- 


,  -  8.589  q12)  (5.32) 

qlf  we  get 

(5.33) 

(5.34) 

(5.35) 

(5.36) 

from  equations 

,  we  get 

-  (5.37) 

-8.859  q12) 

3.0  q2) 


. 


,p  .  d8fcdO.O  - 
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5-3  Computing  Technique 

The  optimization  interval  is  divided  into  24 

sub-intervals,  each  of  one  hour  duration.  The 

required  generation  schedule  is  obtained  by  solving 

equations  (5.15),  (5.26)  and  (5.33).  For  known  values  of  X, 

y  and  Y/oxrw  these  equations  can  be  solved  for  P  , 

(1)0  (2)0  n  T 

q^  and  q^.  By  substituting  for  q^  and  q^  in  equations 
(5.27)  and  (5.16)  the  corresponding  values  of  PH^  and 
P 

K2  can  be  obtained.  The  values  of  A  and  and  Y^)0 

are  found  by  using  the  constraints  of  the  problem,  as 
explained  in  the  following  paragraphs. 

A  appears  in  the  three  scheduling  equations  due 
to  the  point  constraint  that  the  load  demand  is  to  be 
met  at  each  hourly  interval.  That  is,  we  require 


5D  PH1  +  PH2  +  PT 


(5.38) 


This  relationship  must  be  satisfied  at  all  times. 

The  values  of  Y^q  and  Y(2)o  are  chosen  such  that 
the  total  amount  of  water  used  by  each  plant  over  the 
optimization  interval  is  equal  to  the  limited  specified 
quantity.  That  is, 

■T 


q  dt  =  Q. 
0  J  3 


(5.39) 


j  =  1,  2 


. 


Ksrfi  riz.  a  ri  a c  5  s 

sri  1.SV0  :+aslq  rides  bs  iu 
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The  purpose  of  writing  the  computer  program 
is  not  only  to  solve  the  three  scheduling  equations, 
but  also  to  compute  the  values  for  X  and  y  and  Y(2)0* 

The  general  layout  of  the  program  is  given  by  the  flow¬ 
chart  given  on  page  62  of  this  thesis. 

The  scheduling  equations  (5.15),  (5.26)  and  (5.33) 

can  be  solved  for  known  values  of  X,  Y^-jjq  and 
The  constraint  given  by  equation  (5.38)  must  be  satisfied 
at  all  instants  during  the  optimization  interval.  For  this 
reason  a  X  loop  appears  within  the  y  loop  of  the  program. 
The  values  of  y^-jjq  and  Y^)o  can  correcte<^  only  when  a 
trial  schedule  for  24  hours  has  been  calculated. 

Hence,  to  start  with,  some  suitable  values  for  Y^jq 
and  Y(2)o  are  assume(^  and  the  scheduling  equations  are 
subsequently  solved  for  the  24  intervals.  At  each  of  these 
intervals  new  value  value  for  X  is  to  be  obtained.  This 
is  done,  by  assuming  a  suitable  starting  value  for  X 
and  solving  the  scheduling  equations.  Now,  if  the  con¬ 
straint  given  by  equation  (5.38)  is  not  satisfied,  the 
value  of  X  is  increased  by  a  small  increment  and  the 
above  process  is  repeated  until  the  constraint  is  satis¬ 
fied.  Once  this  is  done,  the  computer  leaves  the  loop 
and  goes  on  to  calculate  new  values  for  y^,  Y2 /  and  h^ 
and  in  order  to  compute  the  generation  schedules  for 
the  next  hour.  Glimn  and  Kirchmayer  [14]  have  suggested 
keeping  (t)  and  Z2  (t)  constant  during  all  of  the 


' 

■ 
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FIGURE  5-2 
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optimization  period.  In  the  program  used  here,  new  values 
of  Z1 (t)  and  (t)  are  calculated  for  each  hour.  The  values 
for  y  ^  and  for  a  particular  interval  are  evaluated  by 
using  equation  (5.35)  and  (5.24)  respectively. 

The  integration  is  done  by  simply  adding  the  values 
of  (t)  and  (t)  for  all  the  previous  hours.  Also,  it  is 
assumed  that  during  any  sub-interval  y^  and  y 2  remain  constant. 

,  q£ ,  h^  and  are  also  assumed  to  remain  constant  during 
the  sub-intervals.  For  better  accuracy,  it  may  be  desirable 
to  divide  the  optimization  period  into  a  larger  number  of 
sub-intervals.  Thus,  by  following  the  above-mentioned  proc¬ 
edure,  the  generation  schedule  is  prepared  for  all  the  24 
hourly  sub-intervals. 

To  arrive  at  the  correct  values  for  y^g  and  Y(2)0 
a  method  similar  to  the  relaxation  technique  is  developed 
in  this  thesis.  The  technique  is  a  modification  of  the 
method  suggested  by  Dandeno  [22], 

Suitable  starting  values  of  y^g  an<^  y(2)0 
be  available,  computed  or  guessed  [22],  By  making  rough 
calculations,  based  on  logical  assumptions,  it  is  possible 
to  find  rough  values  of  X,  Y^)o  an<^  Y(2)0‘  Once  suitable 
values  of  y^jq  and  y  (2)  0  are  °ktaine<3'  a  trial  schedule  is 
prepared  by  solving  the  equations  at  each  hourly  interval. 


Y  .’xoi  &&uL&\ 


:•  :  -I 
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Let  the  amount  of  water  used  by  the  two  plants,  over 
the  optimization  interval,  be  and  respectively.  The 
values  of  and  Y(2)0  are  to  correcte<^'  so  that 

the  total  amount  of  water  used  is  and  .  To  do  this 
two  more  schedules  are  prepared  by  changing  the  values  of 


(y 


(DO 


-  Y(2)q^  and  Y  ( 2 )  o  iR  turn.  Hence,  first  Y^-jjq  is 


changed  by  a  small  increment  Ay^,  while  Y^)o  kept 
constant;  and  then  both  Y^q  and  ^(2)0  are  chan9e<^  ty 
Ay 2  so  that  (y^^q  ”  y(2)0^  remains  constant.  Let  the 
water  recruired  by  the  two  plants  in  the  two  cases  be 


I  I 


and  Q''  ;  and  Q,"  '  and  QI'  '  . 


Let  us  define, 


R1 

— 

Qi 

-  Qi 

(5.40) 

R2 

= 

°2 

-  °2 

(5.41) 

AQ{ 

= 

Qi 

-  Qi' 

(5.42) 

AQ' 

= 

°2 

*  °2' 

(5.43) 

AQ{' 

= 

Qi 

- 

(5.44) 

AQ^' 

= 

Q2 

-  Q2" 

(5.45) 

Now, 

let 

the 

required  correction  in 

the 

values 

and  y 

(2)0 

be 

6 y ^  and  6y2»  The  values 

of 

6  y  ^  and 

6y2  are  obtained  by  solving  the  following  equations: 


* 
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(<5y1  - 

«y2) 

aq' 

aqV 

(5, 

,46) 

1  + 

6y 

2 

=  R1 

Ay1 

Ay2 

(<5y1  - 

6y2> 

aq' 

T— -  + 

Ay 

6y2 

AQy 

Ay  2 

-  R2 

(5, 

,47) 

The  corrected  values 

Of 

Y  (1)  0 

and  Y  (2) 0 

are  given  by. 

Y(l)0 

(corrected)  = 

Y 

(1)0  + 

6y1 

(5. 

.48) 

Y  (2)0 

(corrected)  = 

Y 

(2)0  + 

6y2 

(5. 

.49) 

A  new  schedule  is  now  prepared  by  using  the  new 
values  for  Yq^q  and  T^e  Process  maY  have  to  t>e 

repeated  a  number  of  times,  till  the  constraint  given  by 
equation  (5.39)  is  satisfied.  The  number  of  such  iterations, 
which  is  required,  depends  upon  the  extent  to  which  the 
constraint  (given  by  equation  (5 .  39 ))  has  to  be  satisfied 
and  the  starting  values  of  Y^)o  an<^  Y(2)0‘  procedure 

is  found  to  work  quite  satisfactorily  in  all  the  cases 
considered  in  the  following  paragraphs. 

The  complete  computer  program  is  given  in  Appendix 
II.  The  computer  results  for  the  problem  are  given  in 
Appendix  III.  Table  I  on  page  6 6  gives  the  initial  and  the 
corrected  values  for  y^q  and  Y^2)0  and  the  corresPondinP 
total  discharges  of  the  two  hydro-plants. 
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TABLE  I 


Calculation  of  the  Value  of  Y(l)0  and  Y  ( 2 ) 0 


Y  (1)0 

Y  (2)0 

Total  Water  Used 
by  Hydro-Plant 
No.  1  in  24  hrs. 

Total  Water  Used 
by  Hydro-Plant 

No.  2  in  24  hrs. 

Initial  Values 

225 . 000 

65.000 

139.76 

131.17 

1st  Corrected 

241.651 

67.995 

123.62 

153.18 

Values 

2nd  Corrected 

241.259 

68.048 

124.98 

150.05 

Values 

3rd  Corrected 

241.258 

68.050 

125.00 

150.00 

Values 

Amount  of  water  available  for  use  by  plant  No.  1  =  125  K.S.F.-hr. 
Amount  of  water  available  for  use  by  plant  No.  2  =  150  K.S.F.-hr. 
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5 . 4  Solution  by  Alternative  Methods  and  Comparison 
of  Results 

To  test  the  effectiveness  of  the  method  developed 
in  this  thesis,  the  problem  discussed  in  section  5-1  is 
solved  by  using  the  following  methods  as  well. 

Alternative  Method  No.  1: 

Constant  y  Method: 

This  method  has  been  developed  by  Kirchmayer  [18] . 
It  does  not  take  into  account  the  effect  of  variation  of 
head  upon  hydro-plant  characteristics.  Neglecting  trans:- 
mission  loss  terms,  the  scheduling  equations  become, 
for  the  thermal  plant; 


ac 


T 


a  p 


T 


=  X 


(5.50) 


and  for  the  hydro-plant  number  1 


(Yi  -  y2> 


3q. 


ap 


H2 


and  for  hydro-plant  number  2 


aq. 


Y' 


=  X 


ap 


H2 


(5.51) 


(5.52) 


Both  y1  and  y2  remain  constant  during  the  optimi¬ 
zation  interval.  The  problem  is  solved  under  the  same 
assumptions  and  on  the  same  lines  as  those  used  in  section 


5-3. 


■ 


1  ’iTi/fc  Jns^anoo  .ai  s msi  bas 

: 
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Alternative  Method  No.  2: 

One  possible  method  for  operating  the  hydro¬ 
plants  is  to  keep  their  discharge  constant  at  an  average 
value.  The  total  output  of  the  two  hydro-plants  is 
computed  and  the  deficit  is  generated  by  the  thermal 
plant.  The  cost  of  operating  the  system  in  this  way  is 
calculated . 


Alternative  Method  No.  3: 


The  problem  stated  in  section  3-1  is  modified 
by  assuming  that  the  two  hydro-plants  are  located  on 
separate  streams.  The  scheduling  equations  for  the 
thermal  plant  and  the  downstream  hydro-plant  are  the 
same  as  equations  (5.1)  and  (5.2)  respectively.  The 
scheduling  equation  for  the  upstream  plant  becomes: 
(from  equation  (3.47)) 


9P 


HI 


3  P 


Y 


0  9S1 


/ 


HI 


3s. 


dt 


3  P 


(DO 


H  1 


3s. 


(5.53) 


The  three  equations  are  solved  under  the  same 
assumptions  and  by  using  the  same  method  for  numerical 
solution  as  those  discussed  in  section  5-3. 


. 
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Comparison  of  Results 

The  cost  of  operating  the  system  according  to  the 
generation  schedule  obtained  by  using  the  method  developed 
in  this  thesis;  and  by  using  the  alternative  methods  des¬ 
cribed  in  this  section  are  given  in  Table  II  on  page  70 
of  this  thesis.  The  total  amount  of  energy  generated  by 
the  thermal  plant  of  the  system,  over  the  optimization 
interval,  is  also  given  in  that  table. 

It  is  found  that  the  cost  of  operating  the  system 
is  a  minimum  when  we  use  the  scheduling  equations  derived 


in  this  thesis. 


bed  '  r 
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TABLE  II 


Comparison  of  Results 


Method  of  Solution 

Cost  of  Operating 
the  System  for 

24  hours  in  $ 

Total  Mw-hr 
Generated  by 
Thermal  Plant 

Over  24  hours 

1.  tJsing  the  scheduling 
equations  developed  in 
this  thesis 

3209 

729.69 

2.  Constant  y  Method 

3214 

731.16 

3.  Discharge  of  the 

two  plant  kept  constant 
at  an  average  value 

3225 

733.14 

4.  The  problem  is  modified 
by  assuming  that  the 
hydro-plants  of  the 
system  are  on  separate 
streams.  The  problem  is 
solved  by  using  the 
scheduling  equations 
developed  in  section 

3377 

765.18 

3-5. 

■ 
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CHAPTER  VI 

CONCLUSIONS  AND  REMARKS 

General  scheduling  equations  have  been  developed 
for  the  thermal  and  variable  head  hydro-electric  plants 
of  an  interconnected  power  system.  The  scheduling 
equations  for  the  hydro-plants  have  been  obtained  by 
using  the  Euler  equation  in  a  "modified"  form,  developed 
in  this  thesis.  The  transmission  line  losses  have  been 
taken  into  account  while  deriving  the  equations. 

The  scheduling  equations  for  the  hydro-plants  on 
separate  streams  are  shown  to  be  equivalent  to  those  by 
Glimn  and  Kirchmayer  [14]. 

While  deriving  the  scheduling  equations  for  the 
common-flow  plants,  (section  4-4)  it  has  been  assumed 
for  simplicity,  that  only  two  hydro-plants  are  located 
on  the  same' stream.  The  discussion,  however,  is  quite 
general  in  character  and  can  be  extended  to  cases  where 
more  than  two  plants  are  on  the  same  stream.  The  time 
(t)  taken  by  water  to  flow  from  the  upstream  plant  to  the 
downstream  plant  has  been  taken  into  consideration  while 
formulating  the  problem  and  deriving  the  scheduling  equations. 
But,  much  remains  to  be  done  before  the  equations  can  be 
successfully  used.  Suitable  numerical  technique  is  to  be 

ds,  (t)  dsj,(t) 

found  for  evaluating  — : — t-t-t-  ,  -  '  etc.  in  order 

dsk(t)  d§k(t) 


' 

I  i 

1 

*noi;  :ut  9  $ni lu^srfoe  arid  p niviisb  i  ns  rnaldcnq  arid  pnidEi/da^od 
.  ..<••  >:  os  arid  sioi ad  nob  ad  od  aniBrnad  r  .  -  ts  %d  ' 

p.  ■  ;.  ■  /  »  ~  •  1 
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9 S^.  (t-x  ) 

to  find  the  value  of - (discussed  in  section  4-5). 

a  s, 
k 

Also  while  expanding  s^(t-T)  by  Taylor's  series,  it  has 
been  assumed  that  t  is  constant.  However,  strictly 
speaking,  T  depends  upon  the  discharge  of  the  upstream 
plant.  This  further  complicates  the  problem  of  evalua¬ 
ting  *  °kvi°usly/  a  careful  analysis  of  this 

9  s. 

k 

aspect  of  the  problem  is  highly  desirable. 

The  scheduling  equations  obtained  in  this  thesis 
have  been  applied  to  a  small  system  consisting  of  one 
thermal  and  two  hydro-electric  power  plants.  The 
two  hydro-plants  are  located  on  the  same  stream,  t 
and  the  transmission  line  losses  are  assumed  to  be 
zero . 


The  same  problem  has  been  solved  under  the  same 
assumptions  and  simplification  by  using  a  few  other 
(approximate)  methods  described  in  section  5-4.  It  is 
found  that  the  most  economical  generation  schedule  is 
obtained  by  using  the  scheduling  equations  developed  in 


this  thesis. 


■ 
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APPENDIX  1 

System  Characteristics  for  the  Problem  Discussed  in 

Chapter  5 

The  system  consists  of  one  thermal  and  two 

hydro-plants.  The  cost  function  for  the  system  and 

the  characteristics  of  the  hydro-plants  are  taken 

from  reference  [15] .  The  expression  for  the  output 

of  the  hydro-plant  no.  2  was  modified  so  that  the 

incremental  water  usage  curve  (  qj  vs.  P„.)  has  a 

spHj 

positive  slope. 

The  subscripts  1  and  2  refer  to  the  upstream  and 
downstream  plant  respectively. 

1.  Cost  function 

CT  =  0.024  PT2  +  4.0  PT  +  1.0 

2.  Plant  discharges 


When  the  hydro-plants  are  on  the  same  river 


and  when  the  hydro-plants  are  on  separate  streams 
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3.  Effective  Head 


The  tailrace  elevation  and  the  head  losses  in  the 

conduits  are  taken  to  be  constant,  so  that  h^  and  h^  can 

be  represented  by  the  following  relationships: 

t 


h1  =  h1(0) 


0 


qx  at 


When  the  hydro-plants  are  on  the  same  stream 

ft 


^2  ^2  ( 0 ) 


(q0  -  q-,  )  dt 


J  0 


A, 


and  when  the  hydro-plants  are  on  separate  streams 

t 

^2 


h2  =  h2(0) 


q0  dt 


0 


A. 


where,  A^  and  A2  are  the  surface  areas  of  the  reservoirs 
of  the  two  hydro-plants.  These  are  assumed  to  be  constant 
within  the  operating  range,  and  their  values  are, 


A1  =  25  x 


3600  x  10 


43560 


=  2310  Acres 


and 


A2  =  20  x 


3600  x  10 


43560 


=  1651  Acres 


7  7 


4.  Power  Generated 


P 


=  0.064  89q1 (hj-20+38 . 107q  -2 . 86  3q12 ) 


Pr2  =  0 . 076  .  5  q 


5.  Minimum  discharges 


6 .  System  Load  Demand 

The  load  demand  is  assumed  to  remain  constant  at 
300  Mw.  throughout  the  optimization  interval  of  24  hours. 

7.  Total  water  available  for  use 


T 


0 


q±  dt 


125  k . s . f .  hr 


rT 


150  k . s . f .  hr 


0 


where 


T 


24  hours 


9^9flW 
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APPENDIX  2 


Computer  Program 
HYDRO-THERMAL  DISPATCH  COMMON  FLOW  PLANTS 

DIMEN  SI0NXG1(30),XG2(30),PD(30),H1(30),H2(30),D1(30),D2(30),Z1(30) 
X,Z2(30),PT(30),PH1(30),PH2(30),P(30)?ER(30),Z11(30),Z01(30),Z02(3Q 
X),D11(30),C1(30),C2(30),Q1(30),Q2(30),G(30),D12(30),D21(30),D22(30 
X) ,GC(30) ,S1( 10) ,S2( 10) 

READ(5,1)SER,H10.H20,XK1,XK2,IL,JK,KL 

1  F0RMAT(5F5. 1,317) 

READ(5, 11 )Q10,Q20,DG1 ,DG2 ,XG11 ,XG21 ,ERQ,K 
11  FORMAT ( 7  F8 .3,112) 

READ ( 5,2 ) (PD ( I ) ,1=1 ,24) 

2  FORMAT ( 12F5 . 1/12F5. 1 ) 

XG20=XG21 

XG10=XGll 

J  =  1 
L=1 

72  WR I TE ( 6 , 100 ) 

100  FORMAT (1H1) 

WRITE ( 6 , 101 ) 

101  FORMAT (////////) 

WRITE (6,71) 

71  FORMAT (3HHHR,4X,1HL,5X,2HG1,7X,2HG2,5X,2HH1,6X,2HH2,5X,2HD1,4X,2HD 
X2,5X,3HPH1,5X,3HPH2,5X,2HPT,4X,2HPD,5X,3HERR,5X,2HZ1,7X,2HZ2,7X,2H 
XQ1 ,6X  ,2HQ2) 

WRITE(6,102) 

102  FORMAT(/) 

CD=0. 

EF=0. 

IA=1 

G1=XG10 

G2=XG20 

XH1=H10 

XH2=H20 

D029LA=IL,JK,KL 

XL=LA 

A=XL/100000. 

XC1=( ( (G1-G2)/(0.06486*A) )-XHl+20. ) 

XD12=(76.214*76.214-34.356*XC1) 

IF(XD12)29,30,30 
30  XD1 1=SQRT ( XD12 ) 

XD1=(76.2 14+XD1 1)/17.178 

XPHl=0.06486*XDl*(XHl+38.107*XDl-2.863*XDl*XDl-20.) 
XD2=(xH2-5.-(G2/(A*0.076)))/3.0 
I F( XD2 ) 36 ,37,37 

36  XD2=0.0 

37  XPH2=0,076*XD2*(XH2-5.-1.5*XD2) 

XPT=(A-4.0)/0.024 

XP=XPH1+XPH2+XPT 

XER=PD( 1)-XP 
AAR=ABS(XER) 

IF(SER-AAR)9,10,10 


Contd. , 
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9  XD1=(76.214-XD11)/17. 178 
XPHl=0.06486*XDl*(XHl+38.107*XDl-2.863*XDl*XDl-20. ) 
XP=XPH1+XPH1+XPT 
XER=PD(1)-XP 
AAR-ABS(XER) 

IF(SER-AAR)29,10,10 
29  CONTINUE 

10  WRITE(6,4)IA,A,Gl,G2,XHl,XH2,XDl,XD2,XPHl,XPH2,XPT,PD(l),XER,CD,EF 
X ,XD1,XD2 

4  FORMAT( I3,1F6.3,4F8.3,2F6.3,3F8.3,1F6. 1 ,3F8.4,2F8.2) 

Q1(1)=XD1 
Q2(2)=XD2 
HI ( 1 ) =XHl 
H2(2)=XH2 
D 1 ( 1 ) =XD1 
D2(2)=XD2 
Z01 ( 1)=0 
Z02 ( 2 ) =0 

GC(1)=0.012*XPT**2+4,*XPT+1 
D07I=2 ,24 

H1(I)=H1(I-1)-D1(I-1VXK1 

H2( I )=H2( I-l)-D2( I-l)/XK2+XDl/XK2 

Z2( I )=D2( 1-1 )/ ( (H2( I-l)-5.-3,*D2( 1-1) )*XK2) 

Z02 ( I )=Z02 ( I -1 )+Z2 ( I ) 

Zll ( I )=H1( 1-1 )-20. +76. 214*Dl(I-l)-8.589*Dl(l-l)**2. +(0.076/0. 06486 
X)*(H2( I-l)-5.-3.*D2( 1-1) ) 

Z1(I)=D1(I-1)/(Z11(I )*XK1) 

Z01(I)=Z01)I-1)-Z1(I) 

ZG2( I )=XG20/EXP(Z02( I ) ) 

ZG 1  ( I )=XG20/EXP(Z02( I ) ) 

D03LL=IL,JK,KL 

XL=LL 

XL=XL/100000. 

PT ( I )=(XL-4. )/0.024 

C1(I)=(( XG1 ( I ) -XG2 ( I ) )/(0.06486*XL) ) - H 1 ( I )+20. 

D12( I )= (76. 2] 4*76. 214-34. 356*C1( I )0 
IF(D12(I))3, 34,34 
34  Dll (I )=SQRT(D12( I ) ) 

Dl( I ) = ( 76 . 214+D11 ( I ) )/17 . 178 

PHI ( I )=0.06486*D1( I )*(H1( I )+38. 107*D1( I )-2 .863*D1( I )*D1 ( I ) -20 . ) 

C2 ( I  )  =  ( (XG2( I )/0.076*XL) ) -H2( I )  +  5. ) 

D22(I )=C2( I )/3.0 
D2 ( I )=D22 ( I )*(-l ) 

I F ( D2 ( I ) )45,46,46 

45  D2( I )=0 

46  PH2( I )=0.076*D2( I )*(H2(I)-5.-l. 5*D2( I ) ) 

P( I )=PHl( I )+PH2( I )+PT( I ) 

ER( I )  =  PD( I ) - P ( I ) 

AER=ABS(ER( I ) ) 

IF(SER-AER)33,19,19 
33  Dl( I )=(76. 214-D11( I ) )/17. 178 

PH 1 ( I ) =0. 06486* D2( I )*(H1( I )+38. 107*D1( I )-2.863*Dl( I)*Dl(I)-20.) 

P( I )=PH1( I )+PH2( I )+PT( I ) 

ER( I )=PD( I ) - P ( I ) 


Contd. 
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AER=ABS(ER( I ) ) 

IF(SER-AER)3, 19,19 
3  CONTINUE 

19  Q1(I)=Q1(I-1)+D1(I) 

Q2 ( I )=Q2 ( I“1 )+D2( I ) 

G(I)=.012*PT(I)**2+4.0*PT(I)+1. 

GC( I )=GC( I“l)+G( I ) 

WRITE(6 s4) I SXLSXG1 ( I ) SXG2( I ) SH1( I ) ,H2( I ) 9D1( I ) ,D2( I )  ,PH1( I ) ,PH2( I ) 
X,PT(I),PD(I),ER(I),Zl(I),Z2(I),Ql(I),Q2(i) 

7  CONTINUE 
WRITE ( 6  s 103 ) 

103  FORMAT (////) 

WRITE ( 6  9  60 ) 

60  F0RMAT(20X  s4HC0ST) 

C0ST=GC(24) 

WRITE(6  s43)C0ST 
43  F0RMAT(20XS1F1093) 

S 1 ( J ) =Q1 ( 24 ) 

S2(J)=Q2(24) 

R2=Q20-S2.(  I ) 

R1=010-S1(I) 

AR1~ABS(R1) 

AR2=ABS(R2) 

IF(ERQ~AR1 )90s91 S91 
91  IF(ERQ“AR2)90S77  977 
90  IF(L-K)76#76»77 

76  IF(J-2)73,74S75 

73  XG10=XG1 1+DG1 
XG20=XG21 
J=J+1 

GO  TO  72 

74  XG10=XG11+DG2 
XG20=XG21+DG2 
J=J+1 

GO  TO  72 

75  L=L+1 

X1=(S1 (2 )~S1 ( 1 ) )/DGl 
X2=(S1(3)-S1(1))/DG2 
Y1=(S2(2)“S2(1))/DG1 
Y2=(S2(3)“S2(1) )/DG2 
X= ( R1*Y2-R2*X2 I / (X1*Y2  Y1*X2) 

Y=(R1“X*X1)/X2 
XG20=XG21+Y 
XG10=XGll+XfY 
J  =  1 

XG11=XG10 
XG21=XG20 
GO  TO  72 

77  STOP 
END 
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APPENDIX  3 
COMPUTER  RESULTS 

Explanation  of  the  Symbols  in  Computer  Results\ 


Symbols  in  Computer 
Results 


Symbols  used  in  the  remainder 
of  this  thesis  (for  explanation, 
see  page  viii  and  ix.) 
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